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HOLDER REGULARITY 
FOR HYPOELLIPTIC KINETIC EQUATIONS 
WITH ROUGH DIFFUSION COEFFICIENTS 

FRANgOIS GOLSE AND ALEXIS VASSEUR 


Abstract. This paper is dedicated to the application of the DeGiorgi-Nash- 
Moser regularity theory to the kinetic Fokker-Planck equation. This equation 
is hypoelliptic. It is parabolic only in the velocity variable, while the Liouville 
transport operator has a mixing effect in the position/velocity phase space. 
The mixing effect is incorporated in the classical DeGiorgi method via the 
averaging lemmas. The result can be seen as a Holder regularity version of the 
classical averaging lemmas. 


1. The Fokker-Planck equation 

This paper is dedicated to the application of the DeGiorgi method to hypoelliptic 
equations, with rough coefficients. DeGiorgi introduced his technique [7] in 1957 to 
solve Hilbert’s 19th problem. In this work, he proved the regularity of variational 
solutions to nonlinear elliptic problems. Independently, Nash introduced a similar 
technique [23] in 1958. Subsequently, Moser provided a new formulation of the proof 
in [22]. Those methods are now usually called DeGiorgi-Nash-Moser techniques. 
The method has been extended to degenerate cases, like the p-Laplacian, first in the 
elliptic case by Ladyzhenskaya and Uralt’seva [20]. The degenerate parabolic cases 
were covered later by DiBenedetto [8] (see also DiBenedetto, Gianazza and Vespri 
[11, 9, 10]). More recently, the method has been extended to integral operators, 
such as fractional diffusion, in [4, 3] — see also the work of Kassmann [19] and of 
Kassmann and Felsinger [13]. Further application to fluid mechanics can be found 
in [25, 16, 6]. 

Let A = A(t, X, v) be an MA/(R)-valued measurable map on R x R^ x R^ such 
that 

(1) < A{t,x,v) = A{t,x,v)^ < AI 

for some A > 1. Given T > 0, consider the Fokker-Planck equation with unknown 
/ = fit,x,v) G R 

(2) {dt + v- yx)fit,x,v) = dWy{A{t,x,v)yyf{t,x,v)) +g{t,x,v) 
for x,v G R^ and t > —T, where g = g{t,x,v) is given. 
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Assuming that g £ Lfocil—T, oo); x R^)), it is natural to seek / so that 

(3) / G C'((-T,oo);L2(R^ X R^) and V„/G LL([-r, oo); L2(R^ x R^)), 


in view of the following energy inequality: 

( 4 ) 


/ II •) •)IIl2(RNxR'V)^^'S 

J to 


lll/(^) •) •)llL2(R_iVxRJV) -r ^ 

^ lll/(^0, •) •)IIl2(rNxR«‘) + / llff(®Gi •)IIl2(RJVxR'V)II/(Si G •)IIl2(RJVxR'V)^^'S 

Jto 

< 5ll/(^0, •, •)llL2(RiVxRN) + ill5llL2((-T ,oo)xR"xR'^) 

/•* 

ll/(®) ■; •)IIl2(rJVxR'^)^'® • 


'to 


Applying Gronwall’s inequality shows that leads therefore to the following bound 
on the solution of the Cauchy problem for the Fokker-Planck equation with initial 
data /I , G x R^): 

2 2 

ll/(^) ■) OIIlSj-r^NxR") + ^ J ll^«/(Sj G •)llL2(-R^iVxRN‘)'^S 

< (ll/(^0) •) •)IIl2(R'VxR") + II5|Il2((_7’_^)xR"xR")) 


T+t 


for each r > 0 and each t G (—T, t). This bound involves only the bounds on 
the data fL ^ and q. 

This paper is organized as follows. Section 2 establishes a local bound for 
a certain class of weak solutions of the Fokker-Planck equation. The local Holder 
regularity of these solutions is proved in section 3. As in the application of the 
DeGiorgi method to parabolic equations, these two steps involve rather different 
arguments. The main result in the present paper is Theorem 3.1, at the beginning of 
section 3. Yet, the local L°° bound obtained in section 2 is of independent interest 
and is a important ingredient in the proof of local Holder regularity in section 3. 
For that reason, we have stated this local L°° bound separately as Theorem 2.1 at 
the beginning of section 2. 

The arguments used in this paper follow the general strategy used by DeGiorgi, 
with signihcant differences, due to the hypoelliptic nature of the Fokker-Planck 
equation. Earlier results based on Moser’s method are reported in the literature: 
see [24, 26]. The method used in the present paper is especially adapted to kinetic 
models. 

Shortly after completing our proof of local Holder regularity (Theorem 3.1), we 
learned of an independent approach of this problem by Imbert and Mouhot [18]. 
The main difference between [18] and our own work is that Imbert and Mouhot 
follow Moser’s approach, while we follow DeGiorgi’s argument. 


2. The Local L°° Estimate 

Assume henceforth that T > |. All solutions / of the Fokker-Planck equation 
considered here are assumed to satisfy (3) and are renormalized in the sense that, 
for each y G C'^(R) satisfying x(z) = O(z^) as |z| —>■ oo, one has 

(5) {dt + v = div„(AV„x(/)) - x"if)A : (V^/)®^ + 
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in the sense of distributions on (—T, oo) x x R^. 

Notation: for each r > 0, we set 

Q[r] := (—r, 0) x 5(0, r) x 5(0, r). 

The goal of this section is to prove the following local L°° bound. This is the 
first important step in the DeGiorgi method. 

Theorem 2.1. For each A > 1, each 7 > 0, and each q > 12N + 6, there exists 
K = k[N, M, a, 7 , q] G (0,1) satisfying the following property. 

For each Mjq {TV}-valued measurable map A on R x R'^ X R-^ satisfying (1), 
each g G L'^{Q[^]) such that 

and each f G C'{,((—|, 0); L^(5(0, |)^)), solution of the Fokker-Planck equation 
(5t + n • Va;)/= div^(AV„/) + g on (5[|], 
the following implication is true: 

/ f{t,x,v)\dtdxdv < K / < 5 a.e. on Q[\\ ■ 

JqiI] 

The proof of Theorem 2.1 involves several steps, following more or less closely 
DeGiorgi’s original strategy. We shall insist on those steps which significantly differ 
from DeGiorgi’s classical argument. 

2.1. The Local Energy Inequality. Since the solution / of the Fokker-Planck 
equation considered here is renormalized, for each if G C“((—T, 00 ) x R^ x R^), 
one has 

f ff A : (V„x(/) G + f:x"{f){^vf)’^^)dxdvdT 
J-T J JRn xR'V 

xif)idt+v-V/x)'ipdxdvdT-\- f ff gx!{f)f’dxdvdT. 

/r^xr" J-t J 

Since / G Cb{{-T,oo)-, x R^), one can pick a sequence of smooth test 

functions ip converging to a test function of the form ip{T,x,v) = ls<T<t4>{x,v). 
For each (p G G“(R^ x R^) one finds in this way that 

f ff A: (V,x(/)GV,(/) + 0x"(/)(V„/)®")dxdi;dT 
Js JJr'^xr^ 

( 6 ) =// (pxif)is,x,v)dxdv- (pxif){t,x,v)dxdv 

JJR«XR" JJR«XR" 

+ f ff xif)v ■ ^Xpdxdvdr f f f gx'{f)(pdxdvdT . 

Js JJR"xR« Jt JJR«xR" 

Let xif) = ^if ~ c)+ for some c G R and pick rj G C“(R^). Ghoosing the test 
function (p of the form p{x,v) = ri{x)r]{v)'^, we observe that 

^ : (V„x(/) G VxP + 0x"(/)(V„/)®") 

= ? 7 (a;)A : {2r){v){f - c)+V^/ G V^? 7 (v) -k 77(v)^l/>c(V„/)®^) 

since 

xif) = (/ - c)+ and x"(/) = l/>c ■ 
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Hence 


A : (V,x(f) ® + (/-x"(/)(V„/)®") 

= r^(x)A : {2{r]{v)Vy{f - c)+) (g) ((/ - c)+V^? 7 (w)) + vivfC^vif - c)+)®^) 
= r]{x)A : {Vy{r]{v){f - c)+))®^ - r]{x)A : ((/ - c)+Vr 7 (n))®^ , 

since A = . Inserting this identity in (6) shows that 


i // r]{x)r]{vf{f - c)1{t,x,v)dxdv 

JJR«XR" 

j [ [[ r]{x)\Vy{r){v)if - c)+)\‘^dxdvdT 

AJs JjRiVxR" 

Vix)v{vf{f-c)+{s,x,v)dxdv 

JJR"XR« 

+A [ [[ r]{x){f — c)‘^\Vr]{v)fdxdvdT 

"*" / // {f — c)\v ■ Vrj{x)dxdvdT 

+ dif - c)+'q{x)'q{vf dxdvdr . 

Js JJR^xR'V 


(7) 


Remark. The function i® C^, but only with Lipschitz con¬ 

tinuous derivative. Instead of arguing directly with y as above, one should replace 
X by a smooth approximation Xe and passes to the limit as the small parameter 
e-^ 0. 


2.2. The Dyadic Truncation Procedure. This step closely follows DeGiorgi’s 
classical method. For each integer A: > — I, we define 

Tfc :=-i(I + 2-'=), R^:=l(l + 2-^) 

and we set 

Bk := B(0, Rk), Qk := (Tfc, 0) x Bli= Q[Rk]). 

Pick r]k G C'°°(R'^) such that 0 < 77 ^ < I, satisfying 

? 7 fc = lonHfc, r]k = 0onB^_^, and || Vryfclli” < 2 '^+^ . 


Finally, set 


Cfc := i(I- 2 -'=), and/fe = (/- Cfc)+. 



HOLDER REGULARITY FOR HYPOELLIPTIC KINETIC EQUATIONS 
Write inequality ( 7 ) with ij = rjk and c = Ck , for each s G ( Tk - i , Tk )'. one has 

h ' nk { x )' nk { vf ' fl { t , x , v)dxdv 

+ T / // Vk{x)\^v{'nk{v)fk)?dxdvdT 

A JTk dJRNxR" 

<5// r ] kix )' qkivffkis , x , v)dxdv 

J jR'VxR" 

+A / [[ rik { x ) f ^\ Vr ] k { v)fdxdvdr 

JTk-i d jR^xR" 

+ / [[ ^r]k{v)‘^fkV ■Vr]kix)dxdvdT 

JTk-i J Jr'^ xrJv 

[[ gfkVkix)r,kivfdxdvdT. 

1dJ rNxR" 


fTk-i 


Averaging both sides of the inequality above in s G ( Tk - i , Tk ) shows that 


< 


h T]k{x)r]k{vyfkit,x,v)dxdv 

J jR'VxR" 

+ T / 11 rik{x)\Wy{rik{v)fk)\'^dxdvdT 

A JTk djRiVxR" 

2'" f ff r ] k { x ) T ] k { v )'^ fk { s , x , v)dxdvds 

JTk-i d JR^XR'^ 

+A / ff r ] k { x ) fk \ Vrik { v)\‘^dxdvdT 

JTk-^i d jRJVxRW 

+ / // \rik{vfflv-Vrik{x)dxdvdT 

JTk-i d JrJVxR'^ 

+ / // gfkr]k{x)rik{v)'^dxdvdT. 

JTk-i d jRiVxRW 


Set 


Uk ■■= sup i 


Tfc<t<0 JjRJVxR" 


(a;) 77 fc ivYfk {t, X, v)dxdv 


1 

A 


rO 


Tk J Jr^xR^ 


gk{x)\Vy{rikiv)fk)\‘^dxdvdT. 


By construction 


(8) 


0 <Uk < Uk-i < ... < Ui < Uo . 
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2 '" / [[ r]k{x)T]k{vffUs,x,v)dxdvds 

JTk-i j 

[ If Vkix)f^\Vr]kiv)\‘^dxdvdT 

JTk-i J JR'^xR" 


Tk-i J JR'^xR" 


h'HkivY flv ■ Vrik{x)dxdvdT 


< ( 2 '= + A22'=+4 + 1 . . 2 '=+ 2 ) f fldxdvdt 

J Qk-1 

< 2"'=+'(l + 2A) f f^lf^yodxdvdt, 

Qk-1 


SO that 


(9) C/fc < + 2A) / f^lf^>odxdvdt+ f \g\\fk\lfk>odxdvdt. 

jQk-l •'Qk-1 

2.3. The Nonlinearization Procedure. This step starts as in DeGiorgi’s classi¬ 
cal argument. By Holder’s inequality, for each p > 2, 


[ fk^fk>odxdvdt<([ fldxdvdt] |{/fe > 0 } n Qfe-i 
Qk-1 Qk-1 / 


while, assuming that p > 2 is such that i < i — i, 

p z q 


[ \9\\fk\lfk>odxdvdt <\\g\\Li i f fldxdvdt] \{fk > 0} r\ Qk-i\^ r • 
Now, for each A: > 0, 

{fk > 0} = {fk-1 >Ck- Ck-i} = {fk-i > 2-"-'} 
so that, by Bienayme-Chebyshev’s inequality 

\{fk > 0} n Qk-i\< 22'=+2 f fl_,dxdvdt < 22'=+2|Tfc_i|C/fc_i < 3 • 2'^>^+^Uk-i . 

Qk-1 

Hence 


Uk < 22'=+3(1 -b 2A) ( f fPdxdvdt] (3 • 22'=+i[/fc-i)^”i 

yQk-i J 

+ \\9\\li\[ fj^dxdvdt] (3 • 2^'’’+^{7fe_i)^“p“9 
\-JQk-i J 


< 3(1 -b 2A)2^>^+^ulJ f J 


+3j2‘^'^+^ ■ ulJ ’ 


fldxdvdt 


fldxdvdt 
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If one had an inequality of the form 

(11) ( [ fldxdvdA <CkUk-i, 

yQk-i J 

the right-hand side of the inequality above would be the sum of two powers of Uk-i 
with exponents 

2 111 

1 -1-1->1 and l-b- ->1. 

P 2 p g 

In other words, the bound obtained in the present step would result in a nonlinear 
estimate for the linear Fokker-Planck equation. Obtaining a nonlinear estimate for 
the solution of a linear equation is the key of DeGiorgi’s local L°° bound. The 
wanted inequality will be obtained by a variant of the velocity averaging method, 
to be explained in detail below. 


2.4. A Barrier Function. In fact, the velocity averaging method will not be 
applied to fk itself, but to a barrier function dominating fk- Constructing this 
barrier function is precisely the purpose of the present section. Set Fk{t,x,v) := 
fk(t,x,v)rik{x)rikiv)'^. Then, one has 

{dt + v- Vx)Fk - div„(AV„Ffc) = glf^Ck'nk{x)rik{vY + fk'Hkivf'v ■ Vrjkix) 
-2/fcdiv„(?7fc(x)pfc(v)AV?7fc(v)) - 4pfc(a;)?7fc(u)A : Vyfk G Vr7fc(u) - fik 


where fik is a positive Radon measure because the function z i —{z — Ck)+ is convex. 
Set 


Sk := g'i-f>CkVkix)r]kivf -|- fkVkivfv ■ Vrjkix) 
-2 divv{rikix)rikiv) fkAVrjkiv)) - 2r]kix)rikiv)A : Vyfk G Vrik{v) 
and let Gk be the solution of the initial boundary value problem 

{ (^t “t“ u ■ ^x^Gk div^(AV^G/j,) — Sk on C^k—\ ■> 

Gk(t, x,v) = 0 if |u| = or |x| = Rk-i and v ■ x < 0, 

GkiTk-i,x,v) = 0 . 


Hence 

r (dt+v- Vx){Gk - Fk) - dwy{AVx{Gk - Fk)) = /Xfc > 0 on Qk-i, 
■\ (Gk — Fk){t, x,v) = 0 if |u| = Rk-i or |x| = Rk-i and v ■ x < 0, 

[ {Gk-Fk)iTk-ux,v)=0, 
so that, by the maximum principle, 

(13) 0 < Fk < Gk a.e. on Qk-i ■ 


2.5. Using Velocity Averaging. In DeGiorgi original method, the inequality (11) 
follows from the elliptic regularity estimate in the Sobolev space F[^ implied by the 
energy inequality. Together with Sobolev embedding, this leads to an exponent 
p > 2 in (11). 

In the case of the Fokker-Planck equation considered here, the energy inequality 
(4) gives F{^ regularity in the v variables only, and not in (t, x). A natural idea is to 
use the hypoelliptic nature of the Fokker-Planck equation in order to obtain some 
amount of regularity in {t, x). The lack of regularity of the diffusion coefficients, i.e. 
of the entries of the matrix A forbids using the classical methods in Hormander’s 
theorem [17]. 
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There is another strategy for obtaining regularity in hypoelliptic equations of 
Fokker-Planck type, which is based on the velocity averaging method for kinetic 
equations. Velocity averaging designates a special type of smoothing effect for 
solutions of the free transport equation 

{dt+vV,)f = S 

observed for the hrst time in [1, 15] independently, later improved and generalized 
in [14, 12]. This smoothing effect bears on averages of / in the velocity variable v, 
i.e. on expressions of the form 

/ f{t,x,v)(j){v)dv, 

say for (7“ test functions </>. Of course, no smoothing on / itself can be observed, 
since the transport operator is hyperbolic and propagates the singularities of the 
source term S. However, when S is of the form 

S = (i\Vy{A{t,x,v)Vyf) + g 

where g is a given source term in L^, the smoothing effect of velocity averaging 
can be combined with the regularity in the v variable implied by the energy 
inequality (4), in order to obtain some amount of smoothing on the solution / itself. 
A first observation of this type (at the level of a compactness argument) can be 
found in [21]. More recently, Bouchut has obtained more quantitative results, in 
the form of Sobolev regularity estimates [2]. These estimates are one key ingredient 
in our proof. 

By construction 

supp(Gfc) C Qk-i and supp(5'fc) C Qk-i, 

and one has 


{dt + V ■ Vx)Gk — div^(AVt,Gfc) = on R x x R^ . 


By definition 


Sk = Sk,i + div.j, Sk,2 


with 

Sk,i :=glf>CkVkix)r]kiv)‘^ + fkr]kiv)‘^v ■ Vgkix) 

- 2rik{x)rik{v)A : V^/fc (g) Vrjkiv) , 

Sk,2 ■= - ‘2rik{x)rik{v)fkAVrik{v) ■ 

Moreover 


+ + 2'=+3A||V./fc|U.(Q,_,), 


Writing the energy inequality for (12), we find that 

+ \^vGk?dxdvds 

< \\Sk,l\\L'^{Qk_,_)\\Gk\\L^(Q^_^) + ||-S'fc,2||L2(Q^_^)||V„Gfc|lL2(Q^_^) , 
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SO that 

1 /■° 

i sup ||Gfc(t, •, •)|li 2 (B 2 ) + Y \VyGk\'^dxdvds 

Tk-i<t<o '=-1^ ^Jn-1 

< |7fc-ir^^||-5'fc,i||L2(Qj,_i) sup ||Gfc(t,•,-)||i 2 (B 2 ) 

n-i<t<o 

+ ||>S'fc,2||L2(Q^_j)||V^Gfc||L2(Q^_j) 

< |T'fc-i|||«5'fc,i||i2(Q ) + I: sup ||Gfe(t,•,•)||i 2 (B 2 ) 

Tk-i<t<0 

Hence 

1 /■° 

i sup ||Gfc(t, •, •)|li 2 fB 2 ) + T \VyGk\‘^dxdvds 

Tk-i<t<o '=-1^ 

<2|Tfc_i||| 

In particular 

(14) - 4^fc-illl'^fc.illL2(Qj^_i) + 2A|Tfc_i|||S'fe,2|li2(Q^_^), 

The Fokker-Planck equation in (12) is recast as 

{dt + V ■ ^x)Gk = Sk,i + div„(5'fc^2 + AS/^Gk) , 

and we recall that Gfc, S'fc^i and Sk ,2 are supported in Qk-i- 

By velocity averaging, applying Theorem 1.3 of [2] with p = 2, K = n = l,r = 0 
and m = 1, one finds that 

WDl^^Gkh^ + WDl/^Gkh^ 

< Gw||Gfe|li2 + CAr(l + Rk)^^^\\D,Gk\\%\\\Sk,2\\%^ + 

+Gjv(l + + \\Sk,2\\L,^ + A^/^\\DMk\\]!2). 

On the other hand 

WDl/^GkU^ < \\Gk\\%^\\D,Gk\\%^. 


2.6. Using the Sobolev Embedding Inequality. With these inequalities, one 
can estimate 

||uy^Gfc|li2 + WDl^GkWL^ + WDl^GkllL^ 

in terms of Uk-i, as follows. Indeed, 

ll/felli 2 (Qj,) ^ 2|Tfc|[/fc < 3C/fc , and ||V^/fc|l| 2 (Q^) < AUk ■ 

Moreover, by definition of fk ■= (/ — Gk)+, one has 

\\fk\\hiQk-i) ^ ll/fc-llli 2 (Q;,_i) 

and 

|V„/fc| = l/>cjV„/| < l/>Cfc„i|V^/| = |V„/fc_i|, 

so that 

WfkWhiQ^^,) < 3[/fc_i, and ||V^/fc|l| 2 (Q^_^) < AUk-i ■ 
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Thus 

\\Sk,i\\L^ < ||3l/>cJ|L^(Q._,) +2"+'(3V3 + 4A3/2)[/1/2 , 

Besides, for 2 < r < g, one has 

\\9^f>Ck\\L-^(Qk-i) ^ llfflU’'(Qfc_i)l{/ > C'fc} n (5fc_i|2“- 

and we recall that 

\{fk > 0} n Qk-i\ < 22'=+2 f fl_^dxdvdt < 2"'=+"|Tfc_i|C/fe_i < 3 • 2^'^+^Uu-i • 

Qk-1 

Therefore, for 2 < r < g, 

||5fe,i|U^ < ||ff||L-(Q._,)(3-2"'=+'C/fe-i)^-^ +2'=+i(3y3 + 4A3/2)[/i/2 , 
so that 

||^fc,i|li^(Q._0 < 6||g||i.(Q^_^)2(2'=+i)(i-f)c/i:J +22'=+4(27+16A3)C/fe_i. 

||^fc,2lli^(Q._0^3A222'=+6[/,_i. 

Inserting these bounds in the energy estimate (14), we find that 
\\Gk\\h^Q,_,)<9\\Sk,i\\ 

LHQk-i) + 3A||5'fc,2||i2(Q^_^)) 

< 54||g||i.(Q,_,)2(2'=+i)(i-7)[/^^:J + 9(27 + 24A3)22'=+4c/fe_i, 

while 

||Vt,Gfc||^2(Qj^_j) < iA||5'fe,i||^2(Qj^_^) + A^||-S'fe,2|li2(Qj^_^)) 

< 9A||(7llir.(Q,_,)2(2'=+^)(^--)c/fc:J + 3(27A + 24A^)2‘^’^+^Uk-i ■ 

In the inequalities above, one can use Holder’s inequality to estimate ||<7 ||l>-(q^_i) 
as follows: 

ll5l|L’'(Qfe_i) < \\9\\Li{Qk-i)\Qk-l\'^~'^ < IQ[|]I^^^7- 
Summarizing, we have found that 

||Gfe||^2 + ||Vt,Gfe||^2 + ||S'fe,i|1^2 + IIS's,,2 ||l2 < a‘^2'^^{Uf._{ + Uk-i) 

where a = a[A,7] > 0. With the velocity averaging estimate from the previous 
section, this implies that 

||i4y'Gfc||i2 + WDl^CkWl. + WDl^CkWh < b^2^\ul-_f + Uk-i) 
where b = b[N, A, 7] is given by 

52 :=a^(l + (6 + 5A)GA,). 

By Sobolev’s embedding inequality, one has 

\\Gk\\h<Klb^2^'^{ulZ{+Uk-i), 

where Ks is the Sobolev constant for the embedding 

i7i/3(Rx X R^) C LP(Rx X R^), - = - -^—. 

^ ^ p 2 QN + i 

By (13), this implies that 

liF.lli, < Klb‘^2^\ul-J + Uk-i) , i = i - . 
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Using (8), we further simplify the inequality above to obtain 


\Fk\\h < K'U^ + Ul)b^2^'^ulZ{ , 


1 _ 1 1 
p~2~m 


for each k > 1. 


2.7. The Induction Argument. This last step closely follows DeGiorgi’s classical 
argument. We return to (10), and observe that 

/ fj^dxdvdt < / fj!_^dxdvdt = / Fj^_j^dxdvdt 
JQk-l JQk-l •'Qk-l 

<[ FP_^dxdvdt = \\Fk.Z\% ■ 

Qh-o. 


Therefore 


Uk < 3X|(1 + 2A)(1 + Ul)b^ ■ 2^’^+^ulJulJ 

+ 3Ksil + u!)h\\gU, ■ 2^^ul~_t~^ut} 

< 3A§(1 + 2A)(1 + ul)b'^ ■ 2®'=+2Ufel|“" 


F3Ks{l + UZ)n\\g\\L.-2^>^U^4 
where the second inequality above follows from (8). 
With p> 2 being the Sobolev exponent given by 

1 _ 1 1 
p~2~ 6A + 3 ’ 
we choose r = q > 12N + 6 so that 


2 p q r 


2 2 

2 -= 2-1 

p r 


2 2 ^ 2 2 
6 A + 3 ~ g “ ^ 6A + 3 ~ g ^ 


1 _ 1 _ 1 _ 1 _ 1 2 
2 p <7 r 6A + 3 

Besides, in view of (9), one has 


Uo < 8(1 + 2A)k + 7 v^|g[|]|"i < 8(1 + 2A) + =: c[A, 7 ]. 


Thus, setting 


we obtain 


. 1 2 ^ 

01 -.— 1 + . -> 1 , 

6 A + 3 q 


Uk<C- , k>2, 


C[N,A,j,q] := 12A:|(l + 2A)cTO(l + ct)52 + 3A:s(l + ct)H 7 . 

With Vk = U 2 k and p = 2^^(1 + C) — notice that p > 1 — we recast this 
inequality as 

k>l. 

Iterating, we find that 
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k-l 


Elementary computations show that 

k + a{k — 1 ) + a^{k — 2) + ... + a 
= k{l + a + ... + — a(l + 2a + ... + {k — 

„fe-i Q,fc _ 


, a* — 1 c? — 1 , 0 ^ — 1 f ka 

= k -:- a- -— = k -:- a 


a — 1 da a — 1 


a — 1 


a — 1 {a — ly 


a{a^ — 1) — k{a — 1) 


< 


(a-1)2 


so that 


since p > 1. Choosing 


Vk < Eo)“ , for each k>l 


Vo <p 

implies that t4 —>■ 0 as fc —>■ +oo. By dominated convergence, this implies that 


/ {f — V^(^t,x,v)dtdxdv = 0 , 
JQll] 


'Qii] 

i.e. that 

f{t,x,v) < i for a.e. (t,x,v) G Q[\] . 

Because of (9), 

Co < 8(1 + 2K)k + 7v^lQ[i]r/" < (8(1 + 2A) + y|Q[|]|i/2)y^^ 
so that one can choose 

/ _ 2q: \ 


(8(1 + 2A)+7|Q[|]|V2)2 


3. The Local H5lder Continuity 

Our main result in this paper, i.e. the local smoothing effect at the level of 
Holder continuity for the Fokker-Planck operator with rough diffusion matrix, is 
the following statement. 

Theorem 3.1. Let A he an MN(Il)-valued measurable function defined a.e. on 
R X X R^ satisfying (1). Let L he an open interval of R and Lt be an open 
subset o/ R^ X R'^. Let f G C{L; L'^{Q)) and g G L°°(/ x Q) satisfy the Fokker- 
Planck eguation (2) on L x LI. Then there exists a > 0 such that, for each compact 
K C I X Q, one has f G C°’'^{K). 

Notice that, by Theorem 2.1, we already know that / G x LI). 

As in the previous section, the proof of this result follows the general strategy 
of DeCiorgi’s original argument, with significant differences. 

3.1. The Isoperimetric Argument. An important step in the proof of regularity 
in DeCiorgi’s method for elliptic equations is based on some kind of isoperimetric 
inequality (see the proof of Lemma II in [7]). This isoperimetric inequality is a 
quantitative variant of the well-known fact that no Fl^ function can have a jump 
discontinuity. More precisely, given an Fl^ function 0 < m < 1 which takes the values 
0 and 1 on sets of positive measure, DeGiorgi’s isoperimetric inequality provides a 
lower bound on the measure of the set defined by the double inequality 0 < m < 1. In 
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the present section, we establish an analogue of DeGiorgi’s isoperimetric inequality 
adapted to the free transport operator. 

Set Q := X 

Lemma 3.2. Let A > 1 and r/ > 0 be given, and let w G (0,1 — 2“^/-^). There 
exists 9 G (0, i) and a > 0 satisfying the following property. 

Let A = A{t,x,v) be an MnCR)- valued measurable map on R x X R^ 
satisfying (1), and let f,g be measurable functions on QU Q[\\ such that 

{dt+v - Va;)/ = div„(AV„/) + g on QU (5[1] , 

together with 

/ < 1 and [(/I < 1 a.e. on QU (5[1] 

and 

|{/<o}ng|>i|Q|. 

Then 

\{f > l-6>}nQ[w/2]| <r], 

or 

\{0<f <l-9}n{QUQ[l])\>a. 

While DeGiorgi’s isoperimetric inequality is based on an explicit computation 
leading to a precise estimate with effective constants, the proof of Lemma 3.2 is 
obtained by a compactness argument, so that the values of 9 and a are not known 
explicitly. 


Proof. If the statement in Lemma 3.2 was wrong, there would exist sequences A„, 
/„ and gn of measurable functions on Q U Q[l\ satisfying (1) and 

(16) + u • Va;)/„ = div„(AV„/„) + 6f„ onQUQ)!], 

together with 

/„ < 1 and \gn\ < 1 a.e. on Q U g[l] , |{/„ < 0} n Q| > i|Q|, 

and yet 

\{fn > 1 - 2-"} n Q[w/2]| > r ), while |{0 < /„ < 1 - 2“"} n (Q U Q[l])| < 2“" . 


We shall see that this leads to a contradiction. 

First, arguing as in (9), wee see that, for each p G (0,1), there exists Cp > 0 
such that 


|V„/„ 


L2((-|p.0)xB(0.p)D - 


n>2. 


By the Banach-Alaoglu theorem, one can assume that 


fn^f and gn^g in L°°{QUQ[1]), 


while 

Ar,V„fn^h in L2((-|p,0) X B{0,pf) 

for each p G (0,1). 

By the variant of hypoelliptic smoothing based on velocity averaging (Theorem 
1.3 in [2]), one has 

/„ / in Lfoc(Q Q[l]) for all 1 < p < oo . 
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Hence 


(17) {dt + V-V^)/= divyh +g 

together with 

/ < 1 and Igl < 1 a.e. on Q U Q[l], 


on Q U Q[l ], 


|{/<0}nQ| > i|Q|, 


and 

|{/ = 1} n(3[w/2]| > g, while |{0 < / < 1} n (QU(3[1])| = 0. 


One has also 


_0)xB(0.p)2 - ^P' 


But 

(/n)+ -t /+ in L\^c{Q 0 Q[l]) for all 1 < p < oo , 
and, since /« < 1 a.e. on QUQ[1], we conclude that /+ is an indicator function as 
it takes the values 0 or 1 a.e. on Q U <5(1]. Besides, for each p G (0,1), 


I|Vi>/+II^2(_|p_0)xB(0,p) 2 ll^•f^°^“•^IL2(-|p.0)xB(0.p)2 -^P- 

Thus, for a.e. {t,x) G (—|p, 0) x B{0,p), the function v i-A f+{t,x,v) is a.e equal 
0 or 1 and v i-A Vvf+(t,x,v) dehnes an element of L'^(B(0, p)). Hence, for a.e. 
(t,x) G (—|yO, 0) X B(0,p), the function v i-A f+{t,x,v) is either a.e. equal to 0 or 
a.e. equal to 1: see Sublemma on p. 8 in [5]. Therefore, /+ is a.e. constant in the 
V variable. 

Likewise, for each p G (0,1), one has 


II^“(/")+IIl2(-|p,0)xB(0,p) 2 ll^•i'"-°^“•^”IL2(-|p,0)xB(0,p)2 - ^P' 

so that 

dl„V„(/„)+-fc in L2((_3^,0) X B{Q,pf) 
for each p G (0,1). On the other hand, 

(18) {dt+v ■Vx)f+<divvk+l onQUQ)!], 


by convexity of z i-A z+. 

Let us prove that k = 0. For each (f) G C'^{Q U Q[l]), multiplying both sides of 
(16) by 4>{fn)+ and integrating in all variables, one finds that 


j (j)An^vfn ■ Vy{fn)+dtdxdv = j \{fn)\{dt + V ■ V^)4>dtdxdv 
- J {fn)+Anyvfn ' Vv4>dtdxdv + J gn{fn)+4‘dtdxdv 

—>■ J + V ■ Voc)4>dtdxdv — J f^h ■ Vy(j)dtdxdv 

+ J gf+(j)dtdxdv. 

On the other hand, multiplying both sides of (17) by /+</> and integrating in all 
variables, one find that 

J (j)h ■ Vyf+dtdxdv = J ^f^{dt+v-'Vx)4>dtdxdv — J f+h-Vycftdtdxdv 


+ / gf+(j)dtdxdv. 




Since /+ is a.e. constant in the variable v, one has 

Vt,/+ = 0 a.e. on Q U (5[1], 

so that 

k = 0 a.e. on Q U Q[l]. 

Eventually, /+ = lp{t,x) for some measurable P C (—1,0) x B{0, 1), with 

(19) {dt+v-V,)f+<l inP'(QUQ[l]), 
with 

(20) |(Px 5(0,1)) ng| < i|Q| and|(Px5(0,l))nQ[w/2]|>r;. 

Since / is independent of v, the inequality (19) holds in 22'((—1,0) x 5(0,1)) for 
each V G 5(0,1). 

Since /+ is an indicator function, for a.e. (to,xo,vo) G (Q U <5(1]), the function 
s I—/(to + s, xo + svo, Vo) has jump discontinuities. Since 

-^/(io + S,Xo + SVo,Vo) = {dtf + Vo ■ Voof){to + S,Xo + SVo,Vo) < 1, 
as 

one has in fact 

+ s,xo + svo,vo) = (dtf + Vo ■ Va;/)(to + s,xo + svo,vo) < 0. 
as 

On the other hand, if 


— I < fo < —1) |a^ol < 1 — w , and — < t < 0 , |x| < , 

then there exists s > 0 and vo G 5(0, 1) such that {t,x,v) = {to + s,xo + svo,vo). 
Indeed, 

, . ^ 1 ^ 1 ^ 1 1 ^11 \x-xo\ ^ |x| + |xo| Iw+l-w 

t-to > -^a;-(-l) = 1-iw , and \v\ = ^^ ^ -j-= 1. 

t — to t — to i — 2 W 

Therefore 

( 21 ) 

lp(t,x)< lp(to,a;o) for a.e. {to,Xo,t,x)G{-^, -1] x 5(0,1 -w) x (--j, 0) x 5(0,^). 
Since (1 — uj)^ > 1, one has 

|(-|,-l]x5(0,l-w)x5(0,l)|>i|g|. 


so that 


|((-|,-l]x5(0,l-o;))\5|>0. 
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Otherwise, 

i|Q|<|(-|,-l]xi?(0,l-a;)xi3(0,l)| 

= |((-|, -1] X i?(0, 1 - w) X B(0, 1)) n (P X P(0, 1))| < IQ n (P X P(0, 1))| , 

which would contradict the first inequality in (20). 

Choosing (to,a;o) € (—|, —1] x P(0,1 — w)) \ P in (21), we conclude that 

lp(t, x) = 0 for a.e. (t, x) G (—^,0) x P(0, j). 

In other words, one has 

|(PxP(o,i))nQ[|]| =0, 

which is the desired contradiction with (20). □ 

3.2. Zooming in the Fokker-Planck Equation. As in the DeGiorgi original 
proof, the local Holder regularity is obtained by controling the oscillation of so¬ 
lutions of the Fokker-Planck equation on a sequence of domains with shrinking 
diameter. This suggests of course using a zooming procedure based on the scaling 
properties of the Fokker-Planck equation. This step follows the classical DeGiorgi 
argument rather closely. 

For each to € R, xq, vq € and e > 0, we define the transformation 7^[to, xo,vo\ 
by the following prescription: 

%[to, Xo, vo]F{s, y, ^) := F{to + e^s, xq + e^y + e‘^svo,vo + eO ■ 

An elementary computation shows that, if 

{dt+v- V^)F = div„(AV„P) -f G 
then f{s,y,^) = %[to,xo,vo]F{s,y,^) satisfies 

(a^ -h ^ • Vy)f = div^iaV^f) + g 

with 

a{s,y,0 ■= Te[to,xo,vo]A{s,y,^) and g{s,y,^) = e'^%[to,xo,vo]G{s,y,^). 

Observe that a satisfies the same assumption as A, with the same constant A > 1. 

Here is a first application of the zooming transformation defined above. With 
oj chosen as in the previous section, i.e. 0 < w < 1 — , set e = a;/3 in the 

zooming transformation defined above, together with to = 0 and xo = vo = 0. 
Assuming that F satisfies 

(dt + v- V,)P = div„(AV„P) + G on (-^^, 0) X P(0, x P(0, f), 
then /(s, y, ^ = TL/aP, 0, 0]P(s, y, ^ satisfies 

(Ss -h ^ • Vy)/ = div 5 (aV^/) + g on Q[|] 

with 

a(s, y, 0 ■= TL/aP, 0, 0]A(s, y, ^) and g{s, y, ^ = ^7I;/3[0,0,0]G(s, y, C). 

By Theorem 2.1, assuming that |g| < ^ a.e. on Q[|], one has the implication 

/ f^dsdyd^ < k[N, A, oo] ^ a.e. on Q [\]. 

Jq[3/2] 
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In terms of F and G, we arrive at the following statement: assuming that |G| < 1 
a.e. on Q[^], 

f F^dtdxdv < A, oo] ^ F < ^ a.e. on ■ 

Jq[uj/2] 

3.3. Reduction of Oscillation. The second key idea in DeGiorgi’s method for 
proving local regularity is the following important observation, which mixes the 
scaling transformation and the isoperimetric argument. 

Lemma 3.3. There exist /3,/r € (0,1) satisfying the following property. For each 
pair f,g of measurable functions defined a.e. on Q\J (5[1] such that 

(dt+v- Va:)/ = div^(AV^/) + g on QU Q[l] 

with 

I/I < 1 and \g\ < /3 on QU Q[T \, 

one has 

OSCq[(^ 3/54] / < tiOSCQ[(^/2] / . 

Proof. Pick (jj £ (0,1 — 2“^/^), and set 

^:=(|)4A'+2«[7V,A,^,oo]. 

Lemma 3.2 provides us with 9 £ (0, and a > 0. Choose then /3 small enough so 


In - > 


rWI+Ml + 2 ,„i, 

a 9 


Since |/| < 1 on Q U Q[l], one can assume without loss of generality that 

|{/<o}ng|>i|Q|. 


|{/<o}ng| < i|g|, 

we shall argue instead with — / and —g instead of / and g respectively.) 

Consider the sequence of functions defined by induction as follows: 

/fc = - 1) + 1, fo = f- 

One easily check by induction that 

fk < fk-i <---</i</o = /<l a.e. on g U g[l], 

and that fk is a solution of the Fokker-Planck equation on g U g[l] with source 
term 

9k ■■= 0~'^g. 

We shall consider only finitely many terms in this sequence, viz. those for which 


In P 

hTe ■ 


(Notice that the third inequality above follows from the constraint on /3 imposed 
at the begining of this proof.) 

First, one has 

{/ < 0} C {/i < 0} < ... < {fk-i < 0} < {fk < 0} , 
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SO that 


|{/fe<o}nQ|>|{/<o}ng|>i| 


On the other hand 


{fk < 0 } = {fk-i < 0 }u {0 < fk -1 <i-e} 


so that the sequence 


satisfies 


rrik ■■= \{fk < 0} n (Q U Q[l])| 


nik = ruk-i + |{0 < fk-i < 1 - 6»} n (Q U Q[l])| 

k 

= Wo + X! < •^'-1 < 1 - ^} o (<5 u (5[i])|. 

1=1 

It is obviously impossible that 

|{0 < fi-i < 1 — 6*} n (Q U (5[1])| > a for each I = 1,... ,k* , 
Indeed, this would imply that 

\\Q\ + k*a < mo + k*a < trik* < |QI + IQ[1]| > 

which is impossible by our choice of k*. 

Notice that, by our choice of /3, one has 

9~^ P < 1, so that \9~^g\ < I a.e. on Q U Q[l] for A: = 0,..., A:* . 

Applying Lemma 3.2 shows that there exists A: G {0,..., fc* — 1} such that 

\{fk > l-6»}nQ[w/2]| <r]. 

Then 

/ ifj,^^)ldtdxdv = / fl If i^>i-edtdxdv 
JqIu,/2] JqIu,/2] 

— / lf. >i- 9 dtdxdv < rj, 

JQ\co/2] 


SO that 


By definition 


so that 


/fc+i < 5 < 1 - 6* a.e. on . 




/ = 1 + - 1) < 1 - 9’^+^ a.e. on . 


In particular 


OSCq[<^3/54] / < (1 - OSCq[<^/2] / < fJ.OSCQ[^/2] f 


0 < (1 - 6'''*+^ < I. 
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3.4. Implications on Holder Continuity. With the estimates gathered above, 
we conclude the proof of local Holder continuity as follows. 

Observe that, for r > 0 and e G (0,1), 

OSCq[^] T^O, 0, 0]/ = OSCq^j,.] / < OSCq[,^] / , 

where 

( 5 e[r] := (—re^, 0 ] x H( 0 ,re^) x i?( 0 ,re). 

Assume that / and g satisfy the assumptions of Lemma 3.3; then 

Il’^ir2/27[0i 0]ffllL°°(QuQ[l]) — 
for each n > 0. Therefore, Lemma 3.3 implies that 

OSCq[^3/54] T^2/27[0,0,0]/ ^ MOSCq[^/2] T^2/27[0, 0, 0]/ 

< MOSCq[^3/54] T^2~/l7l0, 0, 0]/ 

< /r" OSCq[i^3/54] / < 2^” 

pour tout n > 0. 

In particular 


(a y, 0 e Q[f4] ^ l/((f7)'"A i^r^y, - /(«, 0,0)| < 2^" 


for each n > 0. Therefore 


CJ 


6n+3 


,6n+3 


2.273"+i 
I n other words 


< t < 0, |a;|, |u| < 


lf(t,x,v) - /(O, 0,0)1 < 2/x". 


lf(s,y,0 - /(O,0,0)1 < 2exp 


2.273"+i 

ln(£|max(|s|,| 2 /|,|/|)) 


In — 
27 


In /i 


< ^exp 


ln(^max(|s|,|y|,|^|)) 


In 


In fj, 


27 


2 /2 


Vw 


In /i/ In ^ 


max(|s|,|y|,|e|)'"'^/'''^. 


If / and g belong to L°°{Q U Q;!]) and satisfy the Fokker-Planck equation with¬ 
out satisfying the assumptions of Lemma 3.3 on ||/|lioo(QuQ| 4 ]) and ||3llioo(QuQ[i])) 
replacing / and g respectively with fjL and gjL with 

^ = (1 + II/IIl<»(QuQ[1])) + ^II5IIl“(QuQ[1])^ > 


we conclude that 


with 


l/(s,2/,0 -/(0,0,0)| < C'max(|s|,| 2 /|,|/|)" 


(T:=ln^/ln — 


and 


C := 


2 /2 


/i2 yw 


(1 


Il“(QuQ[1])) (^1+ ^llffllL“(QUQ[l ])j ■ 
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Assume finally that F is a solution of the Fokker-Planck equation with source 
term G on some open neighborhood 17 in R x X R^ of the point {to, xo,vo)- 
Assume further that F,G £ F°°(17). Then 7i[to, a;o, uo]F is a solution of the Fokker- 
Planck equation with diffusion matrix 7i [to,xo, vo\A and source term 7i [to, xo, uo]G. 
Arguing as above with / := 7i[to,xo,vo]F, and setting 

s = t — to , ^ = V — Vo , and y = x — xo — svo , 

we conclude that 


\F{t,x,v) - F{to,xo,vo)\ < G{{1 + |uo|)|t - to\ + |u - vo\ + \x- a;o|)'' 
provided that 


0 < to — t < 


OJ 


54(l+|uo|) ’ 


\x-xo\ < 




54(l + |uo|) ’ 




Since to, and no are arbitrary, this proves that F is locally Holder continuous 
with exponent a. 
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